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ABSTRACT 

Context. The edge of the Central Molecular Zone (CMZj is the location where massive dense molecular clouds with large internal 
velocity dispersions transition to the surrounding more quiescent and lower CO emissivity region of the Galaxy. Little is known about 
the ionized gas surrounding the molecular clouds and in the transition region. 

Aims. To determine the properties of the ionized gas at the edge of the CMZ near Sgr E using observations of N + and C + . 

Methods. We observed a small portion of the edge of the CMZ near Sgr E with spectrally resolved [C n] 158/rm and [N n] 205pm fine 
structure lines at six positions with the GREAT instrument on SOFIA and in [C n] using Herschel HIFI on-the-fly strip maps. We use 
the [N ii] spectra along with a radiative transfer model to calculate the electron density of the gas and the [C n] maps to illuminate the 
morphology of the ionized gas and model the column density of CO-dark H 2 . 

Results. We detect two [Cn] and [Nil] velocity components, one along the line of sight to a CO molecular cloud at -207 km s _1 
associated with Sgr E and the other at -174 km s _1 outside the edge of another CO cloud. From the [N n] emission we find that the 
average electron density is in the range of ~ 5 to 21 cm -3 for these features. This electron density is much higher than that of the 
disk’s warm ionized medium, but is consistent with densities determined for bright diffuse H n nebulae. The column density of the 
CO-dark EL layer in the -207 km s -1 cloud is ~ l-2xl0 21 cnT 2 in agreement with theoretical models. The CMZ extends further out 
in Galactic radius by ~ 7 to 14 pc in ionized gas than it does in molecular gas traced by CO. 

Conclusions. The edge of the CMZ likely contains dense hot ionized gas surrounding the neutral molecular material. The high 
fractional abundance of N + and high electron density require an intense EUV field with a photon flux of order 10 6 to 10 7 photons 
cm~ 2 s _1 , and/or efficient proton charge exchange with nitrogen, at temperatures of order 10 4 K, and/or a large flux of X-rays. Sgr E 
is a region of massive star formation as indicated by the presence of numerous compact H n regions. The massive stars are potential 
sources of the EUV radiation that ionize and heat the gas. In addition X-ray sources and the diffuse X-ray emission in the CMZ are 
candidates for ionizing nitrogen. 

Key words. ISM: clouds — ISM: HII regions — Galaxy: center 


1. Introduction 


The Central Molecular Zone (CMZ), the central ~400 pc by ~ 
80 pc of the Galaxy, is a region filled with massive dense Giant 
Molecular Clouds (GMC) having large internal velocity disper- 


sion. Characteristic cloud line widths are > 20 km s \ compared 
to a few km s 1 for GMCs in the disk (e.g. [Morris & Serabyn 

1996fc [Oka et all 1998'; Lem ere et all 20C 
in the CMZ have higher kinetic tern 

m. In addition, GMCs 
peratures.T^ ~ 30 to 

200 K (Ao et al. 20 ll Mills & Morris 

20131), and densities. 

n(HA > 10 4-5 cm -3 (e.g. Jackson et al. 

1996; Oka et al. 1998; 


iDame et akl 120011: iMartin et alJl2004h . than GMCs in the disk. 
The gas distribution within the inner Galaxy and across the edge 
of the CMZ is thought to be due to the gravitational potential and 
dynamics of the gas there (see overview bv lFerriere et alj[2007l) . 
While the atomic and molecular hydrogen distributions in the 
CMZ have been mapped with spectrally resolved H i and CO, 
respectively, less is known about the distribution and properties 
of the ionized gas in this region. 


The electron distribution of the interstellar medium has been 
determined from dispersion measurements of radio e missi on 
from pulsar s and other compact radi o sou rces (IC ordes & Laziol 
I2002L l2003h . lLazio & Cordesl dl998h and IRovI (l2013h find that 
the average electron density in the inner 100 pc is <n(e)>~ 10 
cm~ 3 . The radial distribution has been modeled as a Gaussian 


with scale length of 145 pc (see iFerriere et al.l 120071) . which 
yield s a value at the edge of the CMZ < n(e)>~ 1 - 2 cm~ 3 . 
Lazio & Cordesl ( 1998 ) proposed that the scattering of radio 
waves is not due to a uniform electron distribution but occurs 
in thin ionized layers at the surface of molecular clouds with 
n(e)>10 3 cm~ 3 and J’/,(e) ~ 10 3 K, or at the interface with the 
Hot Ionized Medium (HIM) with n(e)~ 5-50 cm -3 and T/ ; (e) 
~ 10 5 - 10 6 K. Observations of the spectrally resolved emission 
from ions in the CMZ can help distinguish among these different 
models and provide information about the density, filling factor, 
and ionization mechanisms. The objective of this paper is to ex¬ 
plore the nature of the ionized gas at the edge of the CMZ near 
Sgr E using spectrally resolved emission from the fine structure 
lines, [N 11 ] and [Cn], of ionized nitrogen, N + , and ionized car¬ 
bon, C + , respectively. 

The Balloon-borne Infrared Carbon Explorer (BICE) 
mapped [C 11 ] in the inner Galaxy with modest spatial (15') and 
low spectral (—175 km s 1 ) resolution (Nakag awa et al.l fl998). 
BICE detected bright [Cn] emission from the CMZ from 358/3 
< / < 1° but tapering off outside, and confined to \b\ <1° (see 
their Figures 7 and 8). Within the limitations of BICE’s low 
spectral resolution, iNakagawa et al.l ( 1998 1 found that [C 11 ] has 
a large velocity dispersion towards the CMZ, Av>100 km s 1 
(see their Figure 14), similar to the dispersion seen in l2 CO. In 
addition, the Infrared Space Observatory (ISO) observed [C 11 ] 
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and [Nil] along_a number of lines of sight in the CMZ (see 
iRodrifguez-Fernandez et al i 12004 1 with a spectral resolution ~ 
30 km s _1 . However, much better spectral resolution is needed 
to separate individual components and determine gas character¬ 
istics. 

To understand better the distribution of the ionized gas at 
the edge of the CMZ we observed several lines of sight in 
[Nil] and [Cn] in July 2013 using the heterodyne German 
REceiver for Astronomy at Terahertz frequencies (GREAlQ; 
iHevminck et alJ [20 1 2l) onboard th e Stratospheric Ob servatory 
for Infrared Astronomy (SOFIA) (lYoung et al.ll2012l) . In addi¬ 
tion, prior to the GREAT observations, we made small scale 
strip maps in [Cn] to map the overall distribution of ionized 
gas as part of a Hersche 0 dPilbratt et al.l 1201 Oh Open Time 
2 Programme, using the high spec tral resolution Hetero dyne 
Instrument in the Far Infrared (HIFli fde Graauw et al.ll2010h . 

Here we report on the results of these observations, which 
indicate that there are two regions of highly ionized, high tem¬ 
perature, dense gas: one along the line of sight that includes a CO 
molecular cloud associated with Sgr E at V/ Jr ~ -207 km s 1 and 
the other at V/ sr —174 km s _1 outside the edge of a CO molec¬ 
ular cloud located at larger / in the CMZ. The electron densities 
derived from the [N n] emission, n(e) ~ 5 to 25 cm " 3 , are two 
to three orders of magnitude higher than that in the disk’s warm 
ionized medium (WIM), but consistent with those found from 
analysis of [N n] emission in the Carina Nebula, a luminous H n 
region with numerous O stars (iOberst et al.ll201 lT) . This paper is 
organized as follows. In Section 2 we present our observations, 
while Section 3 analyzes the distribution and density of the ion¬ 
ized gas. Section 4 discusses the possible ionization sources, and 
Section 5 summarizes the results. 



Fig. 1 . A portion of a map of CO integrated intensity from - 
180 to -90 km s 1 of the edge of the CMZ from dFukui et alJ 
120061) . and reprinted with permission from AAAS, based on the 
NANTEN CO(l-O) map which has a beam size ~3'. The region 
mapped with HIFI is outlined by the yellow box and the posi¬ 
tions of the [C 11 ] and [N 11 ] observations with GREAT are indi¬ 
cated by crosses. The approximate CO edge of the CMZ and a 
portion of Loop 1 are indicated, as is the location of Sgr E. 


(Fukui et al. 2006 ). This loop , however, is at a Galactic radius, 
Rg ~670 pc ( Torii et aUl201 ()l) . well outside the CMZ. 


2. Observations of [Cn] and [Nu] 

In Figure Q] we show a_portion of the NANTEN 4-m CO map 
from Figure la in [Fukui et al . ( 1200 6). and reprinted with per¬ 
mission from AAAS, that zooms in on the CO edge of the CMZ. 
The yellow box in Figure Q] outlines the region mapped with 
HIFI, and the positions observed in [C 11 ] and [N 11 ] with GREAT 
are indicated by red crosses. [N 11 ] traces just the highly ionized 
gas, while [C 11 ] traces highly ionized gas and the weakly ion¬ 
ized neutral atomic and molecular gas. Located just inside the 
edge of the CMZ is Sgr E, a region with many radio sources lo¬ 
cated around / ~ 358°7 and b ~ 0° ( iGravIl 19941) . which contains 
18 compact (<2') Hu regions and diffuse emission ( iLisz tlll992l) . 
Sgr E is a region of active star formation asso ciated with a CO 
molecular cloud with V) sr ~ -210 km s _1 (lOka et al Jl 1 998l) . The 
radio re combination line survey, H n Region Discovery Survey 
(HRDS: [Anderson et al.ll200.12012!) detected eight of these Hu 
sources towards Sgr E at V/ sr —208 km s -1 , the same velocity 
as the CO cloud. Along the line of sight outside the CMZ there 
are Giant Magnetic Loops, large scale molecular gas loop-like 
featu res, with large line widths and large sca le flows detected in 
CO dFukui et alJl2006t iFuiishita et al.ll2009l) . The edge of Loop 
1 intercepts the plane at l ~358°2 to 358°4, close to the line of 
sight to the edge of the CMZ, and has a Vi sr —135 km s _1 


1 GREAT is a development by the MPI fur Radioastronomie 
and KOSMA/Universitat zu Koln, in cooperation with the MPI fiir 
Sonnensystemforschung and the DLR Institut fiir Planentenforschung. 

2 Herschel is an ESA space observatory with science instruments 

provided by European-led Principal Investigator consortia and with im¬ 
portant participation from NASA. 


2.1. SOFIA [Cu] and [Nu] data 

We observed the lower fine structure transition of ionized ni¬ 
trogen, [Nu] 3 Pi—> 3 Fo at 1461.1319 GHz (205.2 /mi) and, si¬ 
multaneously, the fine structure line of ionized carbon, [Cn] 
2 Py 2 -* 2 Pi /2 at 1900.5369 GHz (157.7 /an), using GREAT. The 
angular resolution is 15" at [Cn] and 20" at [Nu]. To improve 
the signal-to-noise ratio (SNR) the [C 11 ] and [N 11 ] data were 
smoothed to a velocity resolution of 1.20 km s _1 and 1.57 km 
s' 1 per channel, respectively. The spectra were calibrated using 
the procedure outlined in lGuan et al.1 (120121) . 

The GREAT [C 11 ] observations contained atmospheric fea¬ 
tures in the spectral region observed due to the frequency setting 
for large negative velocities of the source, and the [C 11 ] line was 
in the wing of a strong water line. Moreover, the position switch¬ 
ing observing mode made the presence of this feature a signifi¬ 
cant issue that needed to be addressed in the data reduction. The 
standard atmospheric calibration was not sufficient to produce a 
satisfactory fit to the water feature, and we had to us e a s et of 
non-standard parameters in the atmospheric model dGuan et al.l 
l2012j ) to fit the width and the peak of the water line. The main 
change to the normal procedure was to reduce the contribution 
from the Dry atmosphere in steps, allowing the Wet part (water) 
to fit the atmospheric feature. We implemented a set of six mod¬ 
els, and chose the model that best fit the water feature. This fea¬ 
ture was removed and a zero order spectral baseline fit applied 
to the delivered data products. The resulting spectra were, with 
one exception, fairly regular with residual low order polynomial 
variations and only slight offsets from zero; one spectrum had a 
long period, but well ordered, variation in the baseline. In order 
to fit the baselines better than zero order, the delivered data were 
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then reduced further in CLASS0 and baselines corrected with a 
polynomial fitting routine. Comparison of the GREAT and HIFI 
[Cii] data at one nearly coincident position (see below) shows 
excellent agreement, lending confidence to the removal of the 
water line and overall baseline fittings. The intensities were cor¬ 
rected for a main beam efficiency of 0.67. The rms noise for the 
final [C n] products ranges from 0.16 to 0.26 K per 1.20 km s' 1 
channel, while that for [N n] ranges from 0.09 to 0.23 K per 1.57 
km s' 1 channel. 

We observed six lines of sight (LOS) along b- 0° at lon¬ 
gitudes from Z=358145 to 358175. The reference off-position, 
OFF1,Z=358150 and b -0125 contains weak [C ii] and very weak 
[N ii] emission, shown in Figure [2] so that all the on-positions 
had to be corrected for this emission. To correct the observations 
for emission in OFF1 it was observed separately on two nights 
as an on-position using a different off-position, OFF2,Z=358150 
and Z?=0165. The OFF1 lines appear to be singly peaked and a 
Gaussian fit to the [Cn] emission in OFF1 peaks at V/ sr —195 
km s’ 1 and has a full width half maximum (FWHM) ~ 43 km 
s 1 . The narrow [Cn] feature at V ] sr —58 km s 1 (FWHM ~ 9 
km s 1 ) arises from the 3 kpc arm. A Gaussian fit to the broad, 
weak [N ii] line yields a peak V/ ir —205 km s’ 1 and a FWHM 
-65 km s' 1 , close to the characteristics of the [C ii] line consid¬ 
ering the weakness of the line and the low signal-to-noise ratio 
of this detection. No other [N ii] is detected in the band above 
the 3-cr level. 

As a check on the quality of the GREAT data, given that it 
had to be corrected for the atmospheric water line and emission 
in the OFF1 position, we show in Figure [3] a spectrum from the 
HIFI Galactic Observations of Terahertz (GOT C+) [C n] survey 
(lLanger et al J 120 1 Ol iPineda et al.l 120131) at the closest position, 
l =3581696 and b =010, to the GREAT observation at Z =3581700 
and b =010. The offset in longitude of 14" corresponds to 1.2 
HIFI and 0.9 GREAT beam widths, and some portions of the 
main beams are only about an arcsecond apart. The [C n] emis¬ 
sion from HIFI and GREAT have the same distribution in veloc¬ 
ity, both show a strong peak at —210 km s 1 with intensities 
that generally agree within the uncertainties of the baseline fit¬ 
ting, and both show weak emission centered around ~ -175 km 
s 1 that blends into the stronger feature. The integrated intensi¬ 
ties of these two observations over the velocity range -150 km 
s 1 to -250 km s’ 1 differ by only 10%. Furthermore, the GOT 
C+ spectrum used a different off position than the GREAT ob¬ 
servations, lending further support for our procedure to correct 
the GREAT [C ii] and [N ii] for emission in the GREAT OFF1 
position. 

The [C ii] and [N ii] spectra along b- 0° are shown in 
Figure [4] [C ii] emission was detected at all six positions and 
there are two adjacent, and at times blended, features arising 
from the CMZ in the velocity range —150 km s 1 to -250 km 
s 1 . At the two positions where these components are clearly 
separated, the peaks occur at V/ sr ~ -174 and -207 km s' 1 . In 
general the -174 km s' 1 component is weaker than the -207 
km s' 1 component, and both features are broad with a typical 
FWHM of ~ 25 km s' 1 . No emission is detected at V/ sr —135 
km s' 1 , the velocity of Loop 1. 

[N ii] is only detected at four of the six LOS, and is generally 
much weaker than [Cii] except at Z=358145 where the compo¬ 
nent at V; ST ~ -207 km s' 1 is about half as intense as [C n]. The 
four LOS where [N ii] is detected have an rms noise of -0.1 K 
per 1.57 km s' 1 channel whereas the two lines of sight where 
[N ii] is not detected have a larger rms noise, -0.23 K. The [N n] 


3 http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 2. [Cii] and [Nn] emission in the reference OFF1 position 
at Z=358150 and b-(f.25. The data have been smoothed to 2.4 
km s' 1 for [C ii] and 3.1 km s' 1 for [N n] to improve the signal- 
to-noise. The [Cii] emission peaks at V/ ir —195 km s' 1 while 
that from [N n] peaks - -207 km s' 1 , within the fitting errors of 
the weaker [N ii] signal. The narrow [Cn] emission at V; sr ~ -58 
km s' 1 arises from the 3 kpc arm. [N n] is not detected from this 
arm. 


components were fit with a double Gaussian with one peak at ~ 
-215 km s' 1 and a FWHM -35 km s' 1 and the other component 
is at ~ -170 km s' 1 with a FWHM -25 km s' 1 . All the [Cii] 
features have a large SNR while the [N ii] SNR ranges from -4 
to 20 for the four LOS in which the line was detected. 


2.2. Herschel [Cn] map data 


The GREAT observations provide [C ii] at six and [N n] at 
four longitudinal positions in the plane (b = 0°). To get 
insight into the distribution of [Cn] above and below the 
plane we use [Cul maps made with the high spectral resolu- 
tion HIFI (Ide Gra auw et al.| [2010~) instrument onboard Herschel 
dPilbratt et al.l 20101) . We used a number of On-The-Fly (OTF) 
scans to map the 016 x 014 region shown in Figure[T| All HIFI 
OTF scans were made in the LOAD-CHOP mode using a refer¬ 
ence off-source position at Z=358180, Z?= 1190, about 117 away in 
latitude from the top of the map. The map consists of 13 scans 
in Galactic latitude, centered at b = 0° and stepping in longitude 
at 3' (0105) intervals from 358120 to 358180. The OTF /j-scans 
are 24' long, sampled at every 40 ". At 1.9 THz the angular res¬ 
olution of the Herschel telescope is 12", corresponding to a spa¬ 
tial resolution of -0.5 pc at 8.5 kpc, the distance to the Galactic 
Center. However the effects of scanning rate and the 40" sam¬ 
pling of the OTF maps degrades the effective resolution in b to 
80" while it is 12" in /. Thus spatial structures on scale sizes 
smaller than 80" in b will be strongly attenuated. Furthermore 
the total duration of each OTF scan was typically -2000 sec 
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Fig. 4. [Cn] (blue) and [N n] (red) spectra at the six positions observed with GREAT along b- 050. The Galactic longitude and 
latitude are indicated by GXXX.XX+Y.Y. The two features centered at ~ -174 km s 1 and ~ -207 km s '(which are blended at some 
positions) arise from the edge of the CMZ. [C n] is detected at all 6 positions in the CMZ, but [N n] is only detected at 4 positions. 
The two non-detections in [Nn] are associated with greater rms noise. While [N n] is generally much weaker than [C n], the [N n] 
peak temperature at 1=358°. 45 is about half that of [C n]. 


which provides only a small integration time on each spectrum 
(pixel). Thus the rms in the OTF spectra, 0.22 K in T m b in an 
88" beam and 3 km s ' 1 channel, is larger than that in the corre¬ 
sponding GREAT spectra. 

The OTF map data presented here were processed follow¬ 
ing one of the procedures suggested by the Herschel Science 
Center science staff, by running the HIPE pipeline without an 
off-source subtraction to produce level 1 data. Then the heb- 
Correction routine was applied to the level 1 data to remove the 
HEB standi ng waves and afterwards we followed the procedure 
described in IVelusamv & Fangetl (2 014) to produce the level 2 
data and spectral line maps. The fact that hebCorrection sub¬ 
tracts the matching standing wave patterns from a large database 
of spectra eliminates the need for off-source subtraction. Thus 
in our analysis the processed spectral data are free from any off 
source contamination. We use the processed spectral line scan 
map data to make latitude-velocity (b—V) maps of the velocity 
structure of the [Cn] emission as a function of latitude at each 
of the 13 longitudes of the OTF positions. To increase the sen¬ 


sitivity of the b—V maps we use a velocity resolution of 3 km 
s^ 1 . In all the OTF Z?-V maps the velocity feature near -207 km 
s^ 1 is prominent, but that of the weaker -174 km s 1 feature is 
less so, possibly due to the combined effect of low sensitivity 
and beam dilution. However, the bright [C n] features (> IK) are 
well traced in the OTF maps. To bring out the spatial structure of 
the -207 km s _l velocity feature we generate a [Cn] integrated 
intensity map covering -220 km s 1 to -200 km s 1 . This veloc¬ 
ity range facilitates comparison with the Nobayama CO maps 
of lOka et all (11998l) . which are presented in 10 km s 1 intervals. 
To generate the [C n] map we integrate the intensities in the b- 
V maps at each latitude pixel over the range of velocities from 
-220 km s _1 to -200 km s~ 1 around the [Cn] spectral peak. We 
then use these intensities at all longitudes to generate the 2-D 
map shown in Figure [5] Note that while the intensities are sam¬ 
pled contiguously in latitude, we only have data at 13 discrete 
longitudes with a pixel width of 12". To make the data more vis¬ 
ible on the map we broadened their display in longitude to 60". 
There is little or no [C n] for l < 358:30. 
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Fig. 3. GREAT [Cn] spectrum at the position observed with 
GREAT, Z=3581700 and b-(f.O, that is closest to a [C n] spec¬ 
trum from GOT C+. The GOT C+ position / =3581696 and 
b =010 is offset in longitude by 14" from the GREAT point¬ 
ing, which is about the GREAT beam width. The HIFI data have 
been smoothed to the resolution of the GREAT data, 1.2 km s 1 . 
The [C n] GREAT spectrum is similar to that from HIFI show¬ 
ing both the -207 km s 1 and -174 km s" 1 components. The total 
integrated intensities agree within +10%. 


3. Results 

The relationship between the distribution in [Cn] and CO is 
brought out in Figure [5] where we superimpose on the [C n] strip 
maps the CO i ntegra ted intensity from -210 km s 1 to -220 km 
s 1 (lOka et alJll998h . The only caveat is that the region at the 
edge of the CMZ is irregularly mapped in CO in the Nobayama 
maps. However, the AT&T Bell Lab s CO survey mapped the 
CMZ out to / = ±2° (see Figure 1: iMorris & Serabvniil996 ) 
and also shows a sharp edge in CO emission at / ~35815. 
The -207 km s _I cloud is clearly seen in the position-velocity 
maps from the AT&T B ells Labs CO survey (see Figure 4; 
iMorris & Serabvrj 1996 1. The CO emission comes from the 
cloud associated with Sgr E, an active region of star formation 
with bright compact H ii regions. Though our [C n] map is not 
contiguous in longitude (scans are spaced by 3') the brightest 
emission (white in the color wedge) appears to be an arc above 
the plane slightly outside the lowest CO contour. The [C n] arc 
is most likely limb brightened emission. There does not appear 
to be a similar limb brightening for b <0°. Although we do not 
have any [N n] observations along this arc, we conclude from 
the presence of [N n] along the line of sight to the CO cloud at 
b - 0° that the [N n] and [C n] emission arise from an ionized 
region surrounding the cloud, but there is also a likely contribu¬ 
tion to [C n] from a photon dominated region (PDR) sandwiched 
between the highly ionized hydrogen and dense molecular gas 
traced by CO. Also indicated on the map are 15 of the 18 com- 


0.2 


0.1 

<D 

~o 

¥ ° 

D 

_o 

- 0.1 


Fig. 5. Comparison of the [C n] OTF intensity maps integrated 
over the velocity range -220 km s~* to -200 km s 1 shown as 
coulored strips with the CO map from lOka et all (1 998 ) inte¬ 
grated over the velocity range -220 to -210 km s _I . The CO in¬ 
tegrated intensity maps from lOka et all (119981) are shown only 
in 10 km s' 1 steps, but to improve the signal to noise in the 
[Cn] OTF maps we use a range of 20 km s -1 . The CO map 
integrated from -210 to -200 km s _1 shows the same extent in 
b and the same edge at l ~ 35815 as the -220 to -210 km s _I 
map, but extends further to 35910, and would give the same mor¬ 
phological relationship to the [C n] contours. The color bar at 
the top indicates the [Cn] integrated intensity in K km s 1 . The 
[C n] OTF maps sample in b every 0105 in longitude from 35812 
to 35818, so generating a contour map is impractical due to the 
sparse sampling in /. Instead we have indicated by the white line 
the connection of the peak intensity in each strip. The CO con¬ 
tours, indicated in red, are in steps of 15 K km s _1 . The [C n] is 
strongest slightly beyond the edge of the CO emission above the 
plane. The compact Hu sources towards Sgr E dLisztlll992 ) are 
indicated by white and red X symbols, depending on background 
colour. 

pact H ii sources detected by iLiszt) ( 1992 ) towards Sgr E (the 
remaining three are outside the mapped region). In addition to 
these compact H n regions there appears to be low level diffuse 
emission from l ~ 35818 to 35910 between b = ±011 (see Figure 
8 in lEisztl (11992!) ). The compact Hu sources have electron den¬ 
sities of order several hundred cm -3 and require ionizing stars of 
type B0 or brighter (iLiszdl 1992h . 

All six [C ii] and the four [N ii] spectra in Figure [4] show 
emission in the velocity range -140 to -240 km s _1 . The de¬ 
tection of broad emission in [C n] and [N n] indicates the pres¬ 
ence of warm, and perhaps dense, ionized gas at the edge of 
the CMZ. To gain further insight into the nature of the gas and 
its association with different interstellar medium (ISM) compo¬ 
nents we compare [Cn] and [Nn] to the CO(l-O) and Hi spec¬ 
tra along b = 0° taken from the Three -mm Ultimate Mopra 
Milky Way SurvejQ (ThrUMMS) dBarnes et al.ll201 ll) and the 
Australia Telescope Compact A rray H i Galactic Center surve>0 
(iMcClure-Griffiths et al.ll20T2h . respectively. 

In Figure [6] we plot the CO(l-O) spectra at b = 0° along the 
strip observed with GREAT, Z = 358140 to 358175 and b = 0°. 
The corresponding H i spectra are plotted in Figure [7] The CO 
and H i spectra show that there are two different environments 

4 www.astro.ufl.edu/thrumms 

5 www.atnf.csiro.au/research/HI/sgps/GalacticCenter 
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V (km/s) 

Isr 


Fig. 6. CO(l-O) spectra at b = 0° at the longitudes observed in 
[C ii] and [N n]. The data is taken from the ThrUMMS survey 
(uncorrected for main beam efficiency) and the inset is a key to 
the lines of sight observed. Moving towards the CMZ, larger /, 
CO is not detected until / ~ 358(55 (although there is a hint of 
CO at 358550) at the inner edge of the CMZ, and that only the 
-207 km s' 1 component is present. The -174 km s' 1 CO compo¬ 
nent is detected at 358575, the position furthest from the edge of 
the CMZ and it is considerably weaker there than the -207 km 
s 1 component. 


detected in [Cn] and [N n]. The emission at V/ Jr —207 km s -1 
is associated with both CO and H i, while that for -174 km s -1 is 
associated only with H i until I ~ 358575, where there is a hint of 
CO emission. The distributions of the [C n] and [N n] emission 
for the -207 km s~ 1 component are towards the edge of a CO 
molecular cloud associated with Sgr E, while the -174 km s _1 
distribution is along the edge of another CO molecular cloud, as 
can be seen in the Nobayama CO maps in this velocity range 
(lOka et al.lll998l) . 

To evaluate the properties of the ionized gas we derived the 
parameters of the [C n], [N n], and CO from Gaussian fits to each 
line component. In cases where the lines are blended we derive 
their line parameters using a double Gaussian fit over the veloc¬ 
ity range -100 km s -1 to -300 km s' 1 . The intensities are given in 
Table [I] We calculate the corresponding H i intensities by inte¬ 
grating over the velocity ranges of the two [C n] components be¬ 
cause we cannot get unambiguous Gaussian fits over such broad 
and, at low velocities, flat lines. To avoid double counting the 
intensities we chose -190 km s _1 as the dividing velocity of the 
two components, and their intensities are given in Table □ In 
general, there is a good correlation of [N n] intensity with [C n] 
intensity, which suggests that a significant fraction of the [C n] 
emission arises in the [N n] region. We will use the [C n] and 
[N n] integrated intensities to characterize the properties of the 
ionized gas. 



V (km/s) 

Isr 


Fig.7. Hi spectra at b = 0° along the latitudes observed in 
Cn] and [Nil]. The data are taken from the ATCA Hi survey 
McClure-Griffiths et al. [2012b and the inset is a key to the lines 
of sight observed. Hi is detected at all positions at -207 km s 1 
and is strongest close to Sgr E. Hi is also detected at -174 km 
s" 1 but is weaker than at -207 km s' 1 . 


Table 1. Integrated line intensities 


LOS 

I([C ii])“’ B 

I([N n]) c 

I(CO) 

I(Hi) 

V, jr = -207 km s- 1 

358.45+0.0 

63.4 

26.4 

6.7 

1047 

358.55+0.0 

45.4 

15.4 

26.6 

1904 

358.60+0.0 

47.6 

7.1 

60.6 

2189 

358.65+0.0 

38.6 

- 

101.9 

2539 

358.70+0.0 

45.4 

12.4 

98.9 

2923 

358.75+0.0 

57.4 

- 

110.5 

3235 

V, jr = -174 km s' 1 

358.45+0.0 

21.1 

8.2 

- 

730 

358.55+0.0 

12.2 

8.8 

- 

1212 

358.60+0.0 

43.0 

3.7 

- 

1544 

358.65+0.0 

25.5 

- 

- 

1855 

358.70+0.0 

15.3 

5.3 

- 

2177 

358.75+0.0 

21.8 

- 

9.8 

2522 


a) Integrated intensities are in units of K km s -1 . We only report de¬ 
tections with a SNR > 3, see text, b) Typical rms noise in the [Cii] 
integrated intensity is ~1.4 K km s -1 . c) Typical rms noise in the [N n] 
ntegrated intensity is ~1.3 Kkm s -1 . 


3.1. Radiative transfer solutions for [Cn] and [Nn] 

The relationship of the column density of C + and line inten¬ 
sity, I(rCnD, has been ext ensively discussed in the literature 
(e.g. iGoldsmith et al.l 12012b . The [Cn] peak antenna tempera¬ 
tures reported h ere are «0.3 times the kinetic temperature and, 
as discussed by IGoldsmith et al. ( 2012! ). are in the “effectively 
optically thin limit” in which the column density is propor¬ 
tional to the integrated line intensity. Therefore we have from 
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lLanger et al.l ( 2014 !). Equation A.2, the relationship between the 
column density, N( C + ) in cm -2 and integrated line intensity, 
f T(K)dv , in K km s 1 , 

N(C + ) = 2.92 x 10 15 [1 + 0.5e A£/w '(l + — )]I([CII]) cnT 2 (1) 

n 


where AE/k =91.25K is the energy needed to excite the 
2 -f 3/2 level, n is the density of the collision partner of the 
bulk of the gas, and n cr is the corresponding critical den¬ 
sity. The c ollision rate coefficients for e, H, and H 2 are sum¬ 
marized in Goldsmith et all (20121 ) with an update for EE in 
IWiesenfeld & Goldsmith (2014). Electrons will dominate the 
excitation of C + where we detect [Nn]. In the ionized regions 
where [N 11 ] is detected the temperatures are high enough that 
we can neglect the exponential term and 


N(C + ) = 2.92 x 10 15 [1 + 0.5(1 + 


n cr (e) 

n(e) 


)]/([C7Z]) cnT 


( 2 ) 


The relationship of the column density of ionized nitrogen to 
the integrated line intensity, in units of K km s -1 , in the optically 
thin limit, assuming uniform excitation along the line of sight, is 
given by. 


/ he 3 

T ul (K)dv = ———AuifuNi N + ) cm- 2 , (3) 

%nk(v u i) z 

where A„/ is the Einstein A-coefficient, v„/ the transition fre¬ 
quency, f, is the fractional population of the upper state which 
depends on the density of the collision partner, n, I is in units of 
K km s -1 , and N( N + ) is the column density of ionized nitrogen. 
For the 3 Pi —> 3 Po transition at 205 pm this equation reduces to 

I([N II]) = 5.06 x 10 _16 /iA(N + ) (Kkms‘‘) (4) 

where f\ is the fractional population of N + in the 3 1\ state. If we 
assume that the [N 11 ] emission region is uniform we can replace 
A'(N *) by, n(N + )L, where L is the path length in cm. 



Fig. 8. The intensity of [N 11 ] 205 pm line versus 11 (e) for a range 
of ionized nitrogen column densities; the insert gives the key 
for the column density. The curves flatten out as the density ap¬ 
proaches the critical density. 


of the emitting region and the fractional abundance of ionized ni¬ 
trogen, x(N + ) = n(N + )/n,, where n,= ;z(H)+2h(H2)+h(H + ) is the 
total density of hydrogen nuclei including ions. To estimate n(e) 
we replace the density of ionized nitrogen, n( N + ), in Equation^ 
by ;r(N + )« f , and because the [N 11 ] emission arises from highly, 
or perhaps completely, ionized gas we set n, = n(H + ) = n(e). 
For a uniform density region we can take (V(e) = n(e)L, where 
L is the thickness of the emission region, and can now express 
Equation [5] as 


I ([Nil]) = 5.06 x 10- 16 /i»(N + )L (Kkms- 1 ) (5) 

We can solve for f\ in the optically thin limit from the de¬ 
tailed balance equations for the 3 P 2 , 3 Pi, and 3 Po states along 
with the condition /o + ./1 + f 2 =1. In ionized regions electrons 
dominate the excitation of N + . The collisional de-excitation rate 
coefficient by electrons is given by 

c G(u , [) 

Cui = 8.63 x 10 6 T7TJ’ (6) 

guT"' 5 

where g u is the degeneracy of the upper level and G(u,l) 
is the temperature dependent effective coll ision strength (see 
iLennon & Burke! 1994ilHudson & Bellll2004l) . We calculated the 
electro n collisional rate coefficients using the [Hudson & Belli 
( 2004 ) values for the effective collision strengths. The critical 
density for the 1-0 transition, where collisions out of the 3 Pi 
state balance the radiative transition, A 10 , is n cr (e) ~ 44 cm -3 at 
8000K. In Figure [8] we plot the intensity of the 205 pm line as a 
function of the electron density using Equation [4] for a range of 
ionized nitrogen column densities, iV(N + ) = 10 16 to 10 ls cm -2 . 
The intensity scales linearly with n(e) up to ~ 10 - 15 cm -3 be¬ 
yond which it levels off as n(e) approaches the critical density. 

The [N 11 ] emission comes from ionized gas and so can be 
used to derive the electron density directly if we know the size 


I ([Nil]) = J Twdv = 5.06 x 10 ~ 16 x(N + )Ln(e)fi(n(e)), (7) 

where / is in K km s _1 . Rewriting in convenient units we obtain. 


I([NII]) = 0. l56x- 4 (N + )L pc n(e)fi (n(e)), (8) 

where X -4 is the fractional abundance of N + is in units of 10 4 
and L pc is in pc. Thus we can solve for the electron density, n(e), 
as a function of the observed intensity of the 205 fjm line, given 
the size of the emission region and the fractional abundance of 
N + . At low densities, where «(e)< n cr (e), it can be shown that 


n(e) oc 


I (Nil) 

LpcX- 4(N + ) 


cm 


(9) 


where 0.5 < a < 0.7, for n(e)< n cr (e). In the limit n(e)<< n cr (e), 
such as found in the Warm Ionized Medium in the disk, f\ ~ 
n(e)/n cr (e) and a =0.5. Thus, in general, 11 (e) is only moderately 
sensitive to uncertainties in L and jc(N + ). 


3.2. The -207 km s 1 component 

As seen in Figures [4] and 0 there is [C 11 ] emission along the 
line of sight to the CO molecular cloud with V/„■ —207 km 
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s 1 . The peak in [Cn] intensity for this component is at b > 0° 
at the edge of the lowest contour of CO and is thus consistent 
with limb brightened [C n] emission. Therefore the SOFIA [C n] 
and [N n] spectra along b- 0° likely arise from a layer of gas 
whose size is related to the thickness of the [C n] emission be¬ 
yond the lowest CO contour. The thickness of the [C n] emission 
is about 6' which corresponds to ~15 pc at the distance to the 
CMZ. To solve for the electron abundance in the [N n] regions 
we need, in addition to L pc and /([Nn]), x(N*") appropriate to 
the CMZ. To estimate the fraction of N + in the CMZ we use 
the nitrogen abundance Galactic gradient, 0.07 dex kpc ', for 
Re > 3 kpc, derived from Rollestonetakl (20001) which leads 
to x(N + ) oc exp(-R c /6.2) (see iPineda et al.l 20 13t lLanger et al.1 
120141) . However because we do not have information about how 
it scales to smaller Rq we extrapolate this relationship to the 


half way point, 1.5 kpc. Using the local nitrogen abundance frac¬ 
tion x(N + ) =5.51x10 s from iJensen et al.l (2 007 ) yields x(N + ) = 


1.6xl0~ 4 at Rg= 1.5 kpc. As discussed above, the solutions for 
n(e) are not very sensitive to modest uncertainties in x(N + ). The 
solutions for «(e) and N(N') are given in Table|2] 

The electron density for the -207 km s 1 component varies 
from ~9 to 21 cm 3 for the four LOS with [N n] detections and 
has an average value ~14 cm -3 . Thus the gas traced by [Nn] 
is warm dense ionized gas, significantly denser than that de¬ 
rived for the typica l WIM in the disk, which has n(e) of o rder 
a fewxl0~ 2 cm~ 3 dHaffner et al]l2009tIVelusamv et al.l2012h . In 
the next section we will discuss what might generate and sus¬ 
tain such a high density ionized gas towards the cloud associated 
with Sgr E. 


Table 2. Electron density and nitrogen column density 


V,sr a = 

LOS 

-207 

n(e) b 

-207 

A(N + )" 

-207 

A(H + ) 

-174 

n(e) b 

-174 

N( N + ) 

-174 

N(H + ) 

358.45" 

20.7 

1.5el7 

9.3e20 

9.9 

7.5el6 

4.6e20 

358.55 

14.6 

l.lel7 

6.8e20 

10.4 

7.9el6 

4.8e20 

358.60 

9.1 

6.9el6 

4.2e20 

6.3 

4.8el6 

2.9e20 

358.65 

- 

- 

- 

- 

- 

- 

358.70 

12.8 

9.7el6 

5.9e20 

7.7 

5.9el6 

3.6e20 

358.75 

- 

- 

- 

- 

- 

- 


Average 14.3 l.le!7 6.6c20 8.6 6.5el6 4.0e20 

a) In km s' 1 . b) Densities in cm -3 , c) Column densities in cur 2 and 
rounded to one decimal place, d) All LOS are along b =0°. 


Assuming that the nitrogen is nearly fully ionized, the 
column density of ionized carbon in this region is given by 
N(C + ) =x(C + )/x(N + )N(N + ) For x(C + )/x(N + ) we adopt the value 
of 3.2 dJensen et alJl2007t ISofia et akl 12004) and the values for 
N(C + ) are given in Table [3] along with iV(H + ). The C + column 
densities are of order a fewxlO 17 cm 2 , with an average ~4x 10 17 
cm -2 . 

The [Cn] intensity associated with the [Nn] region, 
/h+([Cii]), is calculated using the inverse of Equation (2), from 
the column density N(C + ), 11 (e) derived from [Nn], and assum¬ 
ing T/, = 8000 K. The exact kinetic temperature is not critical 
as the column density is only weakly dependent on temperature 
through the collisional de-excitation rate coefficient, resulting in 
«„.(e) ~ 507V 1 ' 37 cm' 3 , where 7 4 is the kinetic temperature in 
units of 10 4 K. The results are given in Table[3] 

The [Cn] emission from the neutral PDR, 7h,([Cii]), is the 
difference between the observed total [Cn] intensity, 7([Cn]), 
given in Table Q] and the intensity_calculated from the [N n] re¬ 
gion, 7 h +([Cii]) ( Langer et al.l ’ 2014 1. This quantity, 7 h ,([Cii]) 


= 7([Cn]) - 7h+([Cii]) arises primarily from the CO-dark H 2 
gas (also called the CO-dark gas) and is given in Table [3] 
Negative or very small values of 7 h 2 ([C 11 ]) indicate that al¬ 
most all the [C 11 ] emission arises from the dense ionized gas. 
Unfortunately we do not know the kinetic temperature and hy¬ 
drogen density of the CO-dark iV(H 2 ) to calculate accurately its 
measured C + column density. For illustrative purposes we cal¬ 
culate the column density of C + in the H 2 layer, /V H ,(C t '), and 
the corresponding column density of CO-dark H 2 , A'c (1L)= 
AnUC 1 )/x(C + ) (see lLanger et alJl2014h . assuming 100K and 
300 cm' 3 in Equation (1), along with the appropriate n cr ( H 2 ) 
from lWiesenfeld & Goldsmiith] (120141) . These results are listed in 
Table[3]for all the LOS with positive 7h 2 ([C ii]). Higher densities 
will result in lower column densities and vice versa. 


Table 3. C + column densities and intensities 


LOS" 

N(C + f 

7 H +([Cn]) c 

7n,([Cn]) 

N H ,(C + ) 

w + (H 2 y 

-207" 

358.45 

4.9el7 

64.1 

-0.7 

- 

- 

358.55 

3.5el7 

38.5 

7.0 

4.2el7 

8.2e20 

358.60 

2.2el7 

18.3 

29.3 

1.8el8 

3.5e21 

358.65 

- 

- 

- 

- 

- 

358.70 

3.0el7 

31.3 

14.1 

8.6el7 

1.7e21 

358.75 

- 

- 

- 

- 

- 

Average 

3.4el7 

38.1 

16.8" 

1.0el8 

2.0e21 

-174" 

358.45 

2.4el7 

21.1 

0.1 

3.3el5 

6.4el8 

358.55 

2.5el7 

22.6 

-10.4 

- 

- 

358.60 

1.5el7 

9.9 

33.1 

2.0el8 

3.9e21 

358.65 

- 

- 

- 

- 

- 

358.70 

1.8el7 

13.4 

1.3 

8.2el6 

1.6e20 

358.75 

- 

- 

- 

- 

- 

Average 

2.0el7 

16.9 

11.5" 

7.0el7 

1.4e21 


a) All LOS are at b =0°. b) Column density in cm' 2 and rounded to 
one decimal place, c) Intensity of [C 11 ] in the [N 11 ] emission region in 
K km s' 1 , d) V/ Ir in km s' 1 , e) Negative intensities are not included, f) 
Assumes 7\.=100 K and /tOU) = 300 cm' 3 . 


3.3. The -174 km s 1 component 

The -174 km s' 1 [N 11 ] and [C 11 ] emissions are weaker, and in 
some positions much weaker, than those of the -207 km s 1 
component, and are associated with CO only at / =358575 in 
Figure [6] The CO maps in lOka et all (i!998l) show a CO molec¬ 
ular cloud with an edge at this longitude and extending inwards 
to larger l. The -174 km s' 1 emission is also associated with 
Hi emission (Figure [7} that is weaker than that of the -210 km 
s' 1 cloud. We do not detect this weaker [Cn] component at all 
the positions in the HIFI OTF maps (see discussion in Section 
2.2) and so have no estimate of the size of the [C 11 ] and [N 11 ] 
emission region. For comparison we have calculated densities 
and the column densities of this component assuming the same 
size, 15 pc, for the [N 11 ] emission region as that for the -207 km 
s' 1 cloud. These results are given in Tables [2] and [T] The solu¬ 
tions for the -174 km s' 1 component yield slightly lower average 
electron densities and column densities than for the -207 km s 1 
component. 


4. Discussion 

The analysis of [N 11 ] emission at the edge of the CMZ reveals 
a hot dense plasma with an electron abundance ranging from 
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~ 6 to 21 cm -3 . The densities we derive are much larger than 
those characteristic of the disk’s WIM, but are consistent with 
those derived for very bright nebula with numerous luminous O- 
type stars. For example. lOberst et al.l (1201 it) studied the Carina 
Nebula, which has a very large UV flux, using emission from the 
[N ii] 205 /jm and 122 fim lines along with an excitation modefl, 
and find values of n(e) ranging from a few to over 100 cm -3 . In 
the discussion below we will use the results from our analysis 
of [N n] and [C n] emission to characterize some of the possi¬ 
ble energy sources that can maintain such a hot dense ionized 
plasma. 


4.1. Ionization of Nitrogen 


Nitrogen has an ionization potential of 14.534 eV, higher than 
that of hydrogen at 13.598 eV, and in the ISM can be ionized by 
cosmic rays, X-rays, Extreme Ultraviolet (EUV) photons with 
wavelengths shortward of 853.06 A, electron collisional ion¬ 
ization, or via charge exchange with protons (H + ) at energies 
>0.936 eV, the energy difference between the ionization of H 
and N. In a nearly fully ionized region N + is destroyed primar¬ 
ily via electron radiative and dielectronic recombination, and in 
partially ionized regions also via exothermic charge exchange 
with H atoms. We can estimate the required ionization rate to 
sustain a highly ionized N + gas by calculating the electron re¬ 
combination rate for the solutions for n(e) given in Table [2] The 
reco mbination rate coeffic ients are taken from the UMIST data 
base dMcElrov et al. 20131) and 


Badnelll d2006h for radiative re¬ 


combination and Badnel l et alJ ( 2003 ) for dielectronic recombi¬ 
nation. (In molecular hydrogen clouds N + reacts rapidly with 
FU and initiates a chain of reactions producing nitrogen bearing 
molecules, but this channel is not relevant in the WIM, WNM, 
or H ii regions.) Charge exchange with H has a reaction rate co¬ 
efficient at 8000 K ~5.5xl0" 14 cm 3 s" 1 (see Table 2 of lLin et ail 
( 2005 ;)) and when compared with recombination with electrons 
~8xl0" 13 cm 3 s" 1 , we see that electron recombination domi¬ 
nates for n(e)/n( H) > 0.2. To get an estimate of the degree of 
ionization needed to sustain a highly ionized gas we neglect the 
charge exchange of N + with H in the regime where the gas is 
highly ionized and consider only recombination by electrons. 
For large electron densities in the range 5 to 25 cm" 3 such as 
derived here, the recombination rate ranges from 4xl0" 12 to 
2x10"' 1 s" 1 . Sustaining a high ionization fraction (>0.5) re¬ 
quires an ionization rate comparable to or greater than these val¬ 
ues . For typical WIM conditions in the disk, n(e)~ 10 2 cm" 3 , 
the ionization rate required to keep nitrogen highly ionized is 
two to three orders of magnitude lower. In the following we dis¬ 
cuss the likely sources of ionization at the edge of the CMZ. We 
exclude electron collisional ionization because, due to the high 
ionization potential of N, it requires a very hot electron gas with 
kinetic temperatures of order 10 5 K, a condition met only in the 
HIM portion of the CMZ. 


4.1.1. Cosmic ray ionization 

iGoto et a~fl d2014l) have detected H 3 + in the CMZ and use its 
abundance to derive a lower limit on the hydrogen ionization 
rate, dn > lxlO" 15 s , which they suggest could come from 
cosmic rays or X-rays. If the ionization is due to cosmic rays, 
it is about ten to thirty times larger than the co smic ray ion¬ 
ization rate in the disk (llndriolo & McCaflll2012l) . The cosmic 

6 With two transitions the electron density can be derived without 
having to assume a characteristic size for the emission region. 


ray ionization rate of nitrogen (including the secondary photons 
generated by cosmic rays interacting with hydrogen) is a bout 
a factor of two larger than that of hydrogen ( McElrov et al.l 
12013h . -fewxlO" 15 s" 1 , but is much lower than the values 
needed to offset the recombination by electrons. Although re¬ 
gions of intense cosmic ray ionization may exist with ioni zation 
rates —10 14 s" 1 dBavet et al.ll201 lUMeiierink et al.ll201 ll) . rates 
>10" 12 s" 1 would be required for cosmic rays to explain high 
fractional nitrogen ionization at densities n(e) > 1 cm" 3 . Indeed 
iMeiierink et al . (1201 lh consider rates almost that large in their 
models of photon dominated regions (PDRs) but only consider 
high atomic hydrogen densities, n(H) >10 3 cm" 3 such that the 
fractional ionization is small. Their models cannot be applied to 
the ionized gas probed by our [N n] observations^. 


4.1.2. X-ray ionization 


The CMZ has several different sources of X-rays, includ¬ 
ing an accreting black hole at the center, a bright X-ray 
source IE 1740.7-2942 located l ~ 3 5911 and b —Oil (e.g. 
iHeindletak 1 9931 Wang et alJ 120021) . extended diffuse X-ray 
emission (|Ko^amaet ak|[2007l). as well as over 9000 X-ray point 
sources dMuno et al. 2009) some of which lie at the edge of the 
CMZ between l = 35819 and 35910 around b ~ 0° and could pro¬ 
vide the >10" 15 s" 1 ionization rate of hydrogen required to ex¬ 
plain the H 2 + observations of lGoto et al . ( 2014 ). (Unfortunately, 
these Chandra observations of the CMZ survey the Galactic 
Center do not extend into the Sgr E region.) However, unlike 
ionization by cosmic rays, the cross section for ionization of 
heavy atoms by X-rays can be much larger than hydrogen, pri¬ 
marily due to K shell ionization at X-ray energies >0.4 keV, and 
thus X-rays could play an important role in ionizing nitrogen. 
For example, at 1 keV the ionization cross section of nitrogen is 
7xlO" 20 cm" 2 (IVerner et al.ll 19931) over three orders of magni- 
-ogen (G 


tude larger than that for hydrogen (IGoto et al.ll2014tlWilms et al.l 
2000). Therefore, if the X-ray ionization rate needed to explain 
the production of H 3 + is -fewxlO" 15 s , then the correspond¬ 
ing ionization rate for nitrogen will be of order a fewxlO" 12 s" 1 , 
sufficient to explain the observed large fractional ionization of 
nitrogen at n(e) ~ 5 to 25 cm" 3 . The penetration depth of 1 keV 


X-rays is a column density tV(H+2H2) ~ 10“ cm"“, which is 
sufficient to penetrate the diffuse gas and envelopes of clouds 
in the CMZ (IGoto et al.ll2014l) . Whether X-rays are the primary 
ionization source or not depends on the details of the distribu¬ 
tion and luminosity of the X-ray sources, and we do not have 
sufficient information to make this determination. 


4.1.3. Photoionization 

The photoionization cross section cr for ionizing nitrogen ranges 
from abo ut 10~ 17 to 1.5xl0~ 17 cm 2 from the threshold at ~853 A 
to 500 A (iSamson & Angelll990f) . The total photoionization rate 
is given by f l<853 ^ cr(A)(b(A)dA. where (~4t:J) is the incident 
flux integrated in units of photons cm" 2 s" 1 . It has been esti¬ 
mated that to maintain the ionization of hydrogen in the WIM 
one needs an ioniz ing flux of Lyman c o ntinuum photons >10 5 
photons cm" 2 s" 1 (Ha ffner et alJl2009l) . iRevnolds et al.l (119951) 
estimate an incident Lyman continuum flux 4/r.l ~ 2xl0 6 to 9 

7 Note that lMeiierink et all (2011 ) calculate the N + fractional abun¬ 

dance, x(N + )=/i(N + )/n(H+2H 2 ), as <10" 7 * * * * , which is a low ionization 

scenario. They also calculate the line intensities of the [N n] 205 /j m 

and 122 /im lines, but their predictions fall far below the measured in¬ 

tensities given in Table[I] 


9 










































































W. D. Langer et al.: Ionized gas at the edge of the Central Molecular Zone 


xlO 6 photons cm -2 s 1 , sufficient to explain the ionization in the 
WIM. However nitrogen is ionized by EUV photons beyond the 
Lyman limit, for which the flux must be determined from models 
incorporating a distribution of O and B stars and attenuation by 
the surrounding ISM, including considerations of the dumpiness 
of the ISM (cf. lHaffner et alj[2009l) . In general the flux drops 
sharply beyond the Lyman limit, as can be seen in the model of a 
spectrum of a star cluster in Figure 17 of lKaufman et al.1 (i2006l) . 
where the spectral flux (in ergs/s/Hz/sr/cm 2 ) drops by a factor of 
10 at the Lyman edge and by another factor of 10 by ~400 A. We 
can estimate the flux needed to sustain a high nitrogen ionization 
fraction by balancing photoionization with electron recombina¬ 
tion. 


n(N + )/n(N°) 


L<S53 k ^W)dA 

k r (e)n(e ) 


( 10 ) 


where k r (e ) is the electron recombination rate coefficient, «(N°) 
is the neutral nitrogen number density, and (f> the photon flux 
as a function of wavelength. This equation can be simplified by 
assuming that the ionization cross section is roughly constant 
over the energy range of interest (see above), so «(N + )/h(N°) 
~ (<t>cr)/(k r n(e)), where ® is the total photon flux in cm 2 s" 1 . 
Therefore the fraction of ionized nitrogen is 


f(N + ) = 


®cr 

k r n(e) + ®cr 


(ID 


For example, to sustain /(N + ) = 0.5 in dense ionized gas with 
T=8000 K at n(e) = 1, 10, and 20 cm -3 requires an EUV flux 
~6.5xl0 4 ,6.5x10 s , and 1.3xl0 6 photons cm -2 s" 1 , respectively. 
For a nearly fully ionized nitrogen regime the required intensity 
of the EUV field increases sharply. For example, to maintain 
/(N + ) = 0.9 requires an order of magnitude greater flux than 
for /(N + ) = 0.5, which can be found near young mass ive stel- 
lar clu sters dKaufman et al.1 120061) . iRodrfguez-Fernandez et al.l 
(2004) estimate that the PDRs in the CMZ are illuminated by 
a far-UV radiation field about a factor of 10 3 larger than the lo¬ 
cal ISM, and if this increase extends to the EUV it would be 
sufficient to explain the ionization of nitrogen. One constraint 
on the flux distribution near the -207 km s 1 Sgr E feature 
comes from observations of carbon and helium recombination 
lines. IWenger et al.l (12013l) detected carbon recombination lines 
in two of the eight Sgr E H n sources detected in hydrogen re¬ 
combination lines, but did not detect any helium recombination 
lines. This result indicates that the flux at wavelengths short of 
504.26A (energies > 24.59 eV) is small. 


4.1.4. Proton charge exchange ionization 

The last important ionization mechanism is charge exchange 
with energetic protons. In a highly ionized hydrogen gas, n(H + ) 
> 50 n( H), we can neglect the neutralization of N + by charge ex¬ 
change with H atoms, which yields the ratio of ionized to neutral 
nitrogen. 


f(N + ) = n(N + )/n(N°) 


k x {T)n(H + ) 

k r n(e) 


( 12 ) 


where the charge exchange reaction rate coefficient, k x (T) is a 
Maxwellian average of the cross section over the electron veloc¬ 
ity distribution JLin et al.ll2005 ). The value of k x (T) is an expo¬ 
nential function of temperature because only protons with en¬ 
ergies >0.94 eV can undergo charge exchange. To a very good 


approximation, n(e) - n(H + ) and the fractional abundance of ni¬ 
trogen is 


f(N + ) = 


kAT) 

k x (T) + k r (T)’ 


(13) 


which is independent of n(e) and n(H 1 ). However, the fraction 
/(N + ) is very sensitive to the temperature of the plasma because 
of the endothermic nature of the reaction. The charge exchange 
process for N + + H —> N + H + has been measured (Stebbings 
et al. 1960), but only above ~30 eV, and the low e nergy cross 
sections must be calculated theoretically. [Lin et al.l ( 2005 1 have 
made the most recent calculations of the charge transfer reac¬ 
tion with ground state nitrogen, H + + N H + N + , and use 
these to calculate the rate coefficients, k x for forward and reverse 
reactions. The results for the charge exchange reaction rate co¬ 
efficients are given in their Table 2 and can be used to calculate 
the fractional nitrogen ionization as a function of kinetic temper¬ 
ature. 

However, one must be cautious here, because the theo¬ 
retical rate coefficients derived by iLin et al.l (120051) for H + 
+ N — » N + + H are lower th a n those of Steigman^UilJ 
(1197 111 : iButler & Dalgarno ( j_979 >: IKingdon & Ferlandl (119961 
(see Figure 7 in ILin et ak (2005)) below 2xl0 4 K, which is 
the relevant temperature regime for the WIM and H n regions. 
We solved Equation fl3 l f or the rates derived bv lLin et al.l (2005) 
and lKingdon & Feri and (1996) as a function of temperature and 
these are plotted in Figu re [9] Using the reaction rates from 
iKingdon & Feriand! (119961 ) leads to a high fractional abundance 
at temperatures consistent with those measured in the WIM and 
H ii regions, for temperatures 4000 to 10000 K. However the rate 
coefficients from lLin et al] (120051) lead to a much lower ioniza¬ 
tion fraction of nitrogen at these temperatures, and requires tem¬ 
peratures >15,000 K to ionize nitrogen efficiently via H + charge 
exchange. Thus, depending on which calculation of k x is cor¬ 
rect determines whether or not charge exchange with protons 
can sustain a high ionization fraction of N + in the [N ii] region. 


4.2. Conclusions 

From our analysis of [N ii] emission at the edge of the CMZ near 
Sgr E we find evidence for dense, hot, highly ionized gas. The 
association with [C ii] and CO indicates that this dense ionized 
gas surrounds the PDRs associated with two molecular clouds, 
one with V/ iST —207 km s" 1 and the other —174 km s" 1 . The 
electron densities, n(e) ~ 5 to 21 cm -3 , are considerably higher 
than those in the WIM of the disk, but consistent with those ex¬ 
pected for bright diffuse Hu regions, such as the Carina neb¬ 
ula (|Oberst et akl l201 111 . The electron densities are co nsiderably 
lower than those suggested bv lLazio & Cordesld 19981) . n(e)> 10 3 
cm" 3 , if the photoionized layers of clouds are responsible for the 
scattering of radio waves in the CMZ. However, they are consis¬ 
tent with their other suggestion that the scattering could be due to 
the interface with a hot ionized medium with n(e)~ 5-50 cm -3 
and TAe)~ 10 s - 10 6 K. We have no information about the ki¬ 
netic temperature from our observations, and have assumed typi¬ 
cal WIM kinetic temperatures ~ 8000K. Furthermore, our results 
apply to the ionized regions around the molecular cloud and if 
these are indeed the source of the scattering of radio waves, then 
they occupy a small filling factor in the CMZ and may not be 
representative of the electron density in the bulk of the volume. 

There are three ways to sustain such a dense highly ionized 
hydrogen gas with large nitrogen ionization fraction. One is to 
have a strong EUV radiation field with a flux of order 10 6 to 10 7 
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Fig. 9. The fractional nitrogen abundance as a function of kinetic 
temperature due to production by proton charge exchange and 
destruction by electron recombination. Two charge exchange re¬ 
action rate coefficients are plotte d; the solid line is fromlLin et al.1 
(120051) and the dashed line from lKingdon & Ferlandl (11996 ). 

photons cm~ 2 s _1 , another is to have a high enough kinetic tem¬ 
perature for rapid proton charge exchange with nitrogen, and the 
third is if there is a flux of X-rays sufficient to ionize nitrogen in 
the diffuse gas. In the case of photoionization and X-ray ioniza¬ 
tion it is not necessary for the gas to be very hot. For charge ex- 
chan ge, a gas temperature of or der 5000 K suffices for the rates 
from lKingdon & Ferland! d 1 9961) . but a much higher te mperature 
of ord er 15000 K is required if the rates calculated bv lLin et~ail 
(120051) are correct. We are not able to distinguish which one of 
these processes dominates the ionization of the gas observed in 
[N 11 ], or whether all three contribute. The asymmetry in the [C 11 ] 
limb brightening in b >0° versus b < 0° may be indicative of an 
asymmetry in the distribution of ionizing sources. However we 
lack a complete survey of the Sgr E region in X-rays or EUV 
from which to draw any firm conclusions. 

The edge of the CMZ extends a little farther outwards 
(smaller l ) in ionized components than in CO. The ThrUUMS 
CO map, Figure|6] detects CO at -207 km s 1 out to / = 358°50, 
whereas the GREAT [Cn] and [Nn] observations show emis¬ 
sion out to at least 358°45, a difference of 7 pc. Whether the 
CMZ extends any further cannot be determined as the GREAT 
observations did not extend to smaller values of /, and the HIFI 
OTF maps, which extend to 358?20, are less sensitive than the 
GREAT [C 11 ] data and show only a hint of [C n] at -207 km s _1 
extending to 358°35, a distance of 14 pc from the CO edge. 

We have also estimated the column density of CO-dark H 2 
from [Cn] observations and find an average column density 
of order l-2x!0 21 cm -2 . Models of ma ssive molecular clouds 
dWolfire et al.l 1201 Ot iBolatto et al.l 12013h predict such column 
densities, somewhat independent of the mass of the cloud and of 
the intensity of the UV field. Observations of [C 11 ] and CO to- 
wardsa large ensemble of molecular clouds in the Galactic disk 
dLanger et al.ll20 l4h are consistent with this model. Thus clouds 
in the CMZ and disk have similar CO-dark H 2 layers despite 
differences in UV radiation field and density. 


5. Summary 

We have explored the ionized gas at the edge of the CMZ near 
Sgr E, combining deep integrations of spectrally resolved [C 11 ] 
(158/rm) and [Nn] (205/rm) at six lines of sight along b — 0° 
made with the GREAT instrument on SOFIA and HIFI spectrally 
resolved [C 11 ] on-the-fly strip maps in b at thirteen longitudes 
spaced 0°05 apart. We detect two distinct features in the [Cn] 
and [N 11 ] data, one on the line of sight to a CO molecular cloud 
at V/ sr ~ -207 km s 1 , associated with Sgr E, and the other at ~ 
-174 km s~‘ outside the edge of another CO cloud. 

We find that the brightest emission in the HIFI strip maps 
toward the cloud at -207 km s' 1 forms an arc of [Cn] emis¬ 
sion above the plane at the edge of the lowest CO contours. It 
likely represents limb brightened emission in a highly ionized 
layer outside the molecular cloud. We also detected [N 11 ] emis¬ 
sion at a few positions in this edge, and both [C n] and [N 11 ] 
have very broad lines there with AV of order 25 to 35 km s 1 . 
The electron density in this region was derived from [N n] emis¬ 
sion and a radiative transfer model, assuming a characteristic 
size. We find n(e) ~ 9 to 21 cm~ 3 , about three orders of mag¬ 
nitude higher than that characterizing the disk’s warm ionized 
medium (WIM), much smaller than those of compact H 11 re¬ 
gions, but consistent with densities derived from excitation anal¬ 
ysis of [N 11 ] in bright extended H 11 regions, such as the Carina 
nebula (lOberst et al.l 1201 11) . The electron densities determined 
from the [N 11 ] for the ionized envelope of the cloud associated 
with Sgr E are lower than those suggested to explain the scat¬ 
tering of radio waves used to derive the average electron density 
distribution in the CMZ if the scattering comes from the ionized 
envelope (lLazio & Cordesfl998l) . 

The ionization of such a dense region requires a large flux 
of EUV radiation, and/or X-ray flux, and/or high enough tem¬ 
peratures to allow rapid charge exchange with protons via the 
endothermic reaction, H + + N —> H + N + . The presence of abun¬ 
dant N + requires this region to be hot, and the heating could be 
photoelectron emission from the FUV, EUV, and X-ray radia¬ 
tion fields, dissipation of strong magnetic turbulence, or shock 
heating of supersonic turbulence. This region has many compact 
bright Hu sources, indicating that it is an active massive star 
formation region, and such stars have the capacity to ionize and 
heat the gas in this region. The CMZ is also a source of diffuse 
X-rays and discrete X-ray sources that could provide the X-ray 
flux to ionize the nitrogen (as well as the hydrogen), but it de¬ 
pends on details of the distribution and luminosity of the X-rays. 
Finally, we estimate that this CO cloud has a layer of CO-dark 
H 2 with column density ~l-2xl0 21 cm -2 , consistent with theo¬ 
retical models and observations of clouds in the disk. 

We also detect ionized gas at a velocity —174 km s 1 and 
the lines of sight observed in [C 11 ] and [N 11 ] with GREAT have 
H 1 but no CO except except at / = 358°75 which turns out to be 
the edge of a molecular cloud mapped in CO dOka et alii 199 8). 
However this [C 11 ] emission is weak and the number of [C 11 ] 
detections in the HIFI map is insufficient to determine the size 
of the ionized envelope for this cloud. Therefore to estimate the 
properties of the envelope we assumed the same thickness as 
that for the -207 km s -1 cloud. An analysis of the [N 11 ] emission 
using parameters similar to those for the -207 km s 1 cloud sug¬ 
gests that it too is a hot high density ionized gas with n(e)~ 6 to 
10 cm -3 . The physical conditions in this ionized component are 
likely to be similar to those at the edge of the molecular cloud 
at —207 km s 1 and have similar heating and ionizing sources. 
Larger scale maps of spectrally resolved [C 11 ] and [N n] will be 
needed to understand further the role of the ionized gas at the 
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edge of the CMZ. Finally, the CMZ is seen to extend further out 
in ionized gas from the CO edge by upwards of at least 7 to 14 
pc. 
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